Insights from Genomes and Genetic Epidemiology of SARS-CoV-2 isolates from the state of Andhra Pradesh

Coronavirus disease (COVID-19) emerged from a city in China and has now spread as a global pandemic affecting millions of individuals. The causative agent, SARS-CoV-2 is being extensively studied in terms of its genetic epidemiology using genomic approaches. Andhra Pradesh is one of the major states of India with the third-largest number of COVID-19 cases with limited understanding of its genetic epidemiology. In this study, we have sequenced 293 SARS-CoV-2 genome isolates from Andhra Pradesh with a mean coverage of 13,324X. We identified 564 high-quality SARS-CoV-2 variants, out of which 15 are novel. A total of 18 variants mapped to RT-PCR primer/probe sites, and 4 variants are known to be associated with an increase in infectivity. Phylogenetic analysis of the genomes revealed the circulating SARS-CoV-2 in Andhra Pradesh majorly clustered under the clade A2a (94%), while 6% fall under the I/A3i clade, a clade previously defined to be present in large numbers in India. To the best of our knowledge, this is the most comprehensive genetic epidemiological analysis performed for the state of Andhra Pradesh.

Intra-host variability in global SARS-CoV-2 genomes as signatures of RNA editing: implications in viral and host response outcomes

During the course of the COVID-19 pandemic, large-scale genome sequencing of SARS-CoV-2 has been useful in tracking its spread and in identifying Variants Of Concern (VOC). Besides, viral and host factors could contribute to variability within a host that can be captured in next-generation sequencing reads as intra-host Single Nucleotide Variations (iSNVs). Analysing 1, 347 samples collected till June 2020, we recorded 18, 146 iSNV sites throughout the SARS-CoV-2 genome. Both, mutations in RdRp as well as APOBEC and ADAR mediated RNA editing seem to contribute to the differential prevalence of iSNVs in hosts. Noteworthy, 41% of all unique iSNVs were reported as SNVs by 30th September 2020 in samples submitted to GISAID, which increased to [~]80% by 30th June 2021. Following this, analysis of another set of 1, 798 samples sequenced in India between November 2020 and May 2021 revealed that majority of the Delta (B.1.617.2) and Kappa (B.1.617.1) variations appeared as iSNVs before getting fixed in the population. We also observe hyper-editing events at functionally critical residues in Spike protein that could alter the antigenicity and may contribute to immune escape. Thus, tracking and functional annotation of iSNVs in ongoing genome surveillance programs could be important for early identification of potential variants of concern and actionable interventions. GRAPHICAL ABSTRACT O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=177 SRC="FIGDIR/small/417519v3_ufig1.gif" ALT="Figure 1"> View larger version (41K): org.highwire.dtl.DTLVardef@12b6ac2org.highwire.dtl.DTLVardef@16df897org.highwire.dtl.DTLVardef@dbbec2org.highwire.dtl.DTLVardef@c8de14_HPS_FORMAT_FIGEXP M_FIG C_FIG

Initial insights into the genetic epidemiology of SARS-CoV-2 isolates from Kerala suggest local spread from limited introductions

Coronavirus disease 2019 (COVID-19) rapidly spread from a city in China to almost every country in the world, affecting millions of individuals. Genomic approaches have been extensively used to understand the evolution and epidemiology of SARS-CoV-2 across the world. Kerala is a unique state in India well connected with the rest of the world through a large number of expatriates, trade, and tourism. The first case of COVID-19 in India was reported in Kerala in January 2020, during the initial days of the pandemic. The rapid increase in the COVID-19 cases in the state of Kerala has necessitated the understanding of the genetic epidemiology of circulating virus, evolution, and mutations in SARS-CoV-2. We sequenced a total of 200 samples from patients at a tertiary hospital in Kerala using COVIDSeq protocol at a mean coverage of 7,755X. The analysis identified 166 unique high-quality variants encompassing 4 novel variants and 89 new variants identified for the first time in SARS-CoV-2 samples isolated from India. Phylogenetic and haplotype analysis revealed that the circulating population of the virus was dominated (94.6% of genomes) by three distinct introductions followed by local spread, apart from identifying polytomies suggesting recent outbreaks. The genomes formed a monophyletic distribution exclusively mapping to the A2a clade. Further analysis of the functional variants revealed two variants in the S gene of the virus reportedly associated with increased infectivity and 5 variants that mapped to five primer/probe binding sites that could potentially compromise the efficacy of RT-PCR detection. To the best of our knowledge, this is the first and most comprehensive report of genetic epidemiology and evolution of SARS-CoV-2 isolates from Kerala.

High throughput detection and genetic epidemiology of SARS-CoV-2 using COVIDSeq next generation sequencing

The rapid emergence of coronavirus disease 2019 (COVID-19) as a global pandemic affecting millions of individuals globally has necessitated sensitive and high-throughput approaches for the diagnosis, surveillance and for determining the genetic epidemiology of SARS-CoV-2. In the present study, we used the COVIDSeq protocol, which involves multiplex-PCR, barcoding and sequencing of samples for high-throughput detection and deciphering the genetic epidemiology of SARS-CoV-2. We used the approach on 752 clinical samples in duplicates, amounting to a total of 1536 samples which could be sequenced on a single S4 sequencing flow cell on NovaSeq 6000. Our analysis suggests a high concordance between technical duplicates and a high concordance of detection of SARS-CoV-2 between the COVIDSeq as well as RT-PCR approaches. An in-depth analysis revealed a total of six samples in which COVIDSeq detected SARS-CoV-2 in high confidence which were negative in RT-PCR. Additionally, the assay could detect SARS-CoV-2 in 21 samples and 16 samples which were classified inconclusive and pan-sarbeco positive respectively suggesting that COVIDSeq could be used as a confirmatory test. The sequencing approach also enabled insights into the evolution and genetic epidemiology of the SARS-CoV-2 samples. The samples were classified into a total of 3 clades. This study reports two lineages B.1.112 and B.1.99 for the first time in India. This study also revealed 1,143 unique single nucleotide variants and added a total of 73 novel variants identified for the first time. To the best of our knowledge, this is the first report of the COVIDSeq approach for detection and genetic epidemiology of SARS-CoV-2. Our analysis suggests that COVIDSeq could be a potential high sensitivity assay for detection of SARS-CoV-2, with an additional advantage of enabling genetic epidemiology of SARS-CoV-2.

Analysis of the potential impact of genomic variants in SARS-CoV-2 genomes from India on molecular diagnostic assays

An isolated epidemic of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) causing Coronavirus Diseases (COVID-19) originating in Wuhan, China has now rapidly emerged into a global pandemic affecting millions of people worldwide. Molecular detection of SARS-CoV-2 using reverse transcription polymerase chain reaction (RT-PCR) forms the mainstay in screening, diagnosis and epidemiology of disease. The virus has been evolving through base substitutions. The recent availability of genomes of SARS-CoV-2 isolates from different countries including India motivated us to assess the presence and potential impact of variations in target sites for the oligonucleotide primers and probes used in molecular diagnosis. We catalogued a total of 132 primers or probes sequences from the literature and the public domain. Our analysis revealed a total of 125 unique genetic variants in 80 either primers or probes binding sites. A total of 13 unique variants had allele frequency of [≥] 1% in Indian SARS-CoV-2 genomes mapped to the primers or probes binding sites. A total of 15 primers or probes binding sites had cumulative variant frequency of [≥] 1% in the SARS-CoV-2 genomes. These included primers or probes sites which are widely used in India and across the world for molecular diagnosis as well as approved by national and international agencies. This highlights the need for sequencing genomes of emerging pathogens to make evidence based policies for development and approval of diagnostics. To the best of our knowledge, ours is the most comprehensive analysis of genomic variants in genomes of SARS-CoV-2 isolates from India and their potential impact on efficacy of molecular diagnostics.

Tapestry: A Single-Round Smart Pooling Technique for COVID-19 Testing

The COVID-19 pandemic has strained testing capabilities worldwide. There is an urgent need to find economical and scalable ways to test more people. We present Tapestry, a novel quantitative nonadaptive pooling scheme to test many samples using only a few tests. The underlying molecular diagnostic test is any real-time RT-PCR diagnostic panel approved for the detection of the SARS-CoV-2 virus. In cases where most samples are negative for the virus, Tapestry accurately identifies the status of each individual sample with a single round of testing in fewer tests than simple two-round pooling. We also present a companion Android application BYOM Smart Testing which guides users through the pipetting steps required to perform the combinatorial pooling. The results of the pooled tests can be fed into the application to recover the status and estimated viral load for each individual sample. NOTE: This protocol has been validated with in vitro experiments that used synthetic RNA and DNA fragments and additionally, its expected behavior has been confirmed using computer simulations. Validation with clinical samples is ongoing. We are looking for clinical collaborators with access to patient samples. Please contact the corresponding author if you wish to validate this protocol on clinical samples.